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Abstract 

The Einstein Equivalence Principle has as one of its implications that the 
non-gravitational laws of physics are those of special relativity in any local 
freely-falling frame. We consider possible tests of this hypothesis for systems 
whose energies are due to radiative corrections, ie. which arise purely as a 
consequence of quantum field theoretic loop effects. Specifically, we evaluate the 
Lamb shift transition (as given by the energy splitting between the 251/2 
2 Pi /2 atomic states) within the context of violations of local position invariance 
and local Lorentz invariance, as described by the THe^ formalism. We compute 
the associated red shift and time dilation parameters, and discuss how (high- 
precision) measurements of these quantities could provide new information on 
the validity of the equivalence principle. 
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I Introduction 


The Einstein Eqnivalence Principle (EEP) is fonndational to onr nnderstanding of 
gravity. It states that (i) all test bodies fall with the same acceleration regardless 
of their composition (the weak eqnivalence principle, or WEP) and (ii) the ontcome 
of any local nongravitational test experiment is independent of the velocity and the 
spacetime orientation and location of the (freely-falling) apparatns . Theories which 
obey the EEP, such as general relativity and Brans-Dicke Theory are called metric 
theories because they endow spacetime with a metric that couples universally to 
all non-gravitational helds. Non-metric theories do not have this feature: they break 
universality by coupling auxiliary gravitational fields directly to matter. In this con¬ 
text a violation of the EEP means the breakdown of either Local Position Invariance 
(LPI) or Local Lorentz Invariance (LLI) (or both) so that observers performing local 
experiments could detect effects due to their position (if LPI is violated) or their 
velocity (if LLI is violated) in an external gravitational environment by using clocks 
and rods of differing composition. Limits on LPI and LLI are set by gravitational 
red-shift and atomic physics experiments respectively |^, |^, |^ , each of which compares 
relative frequencies of transitions between particular energy levels that are sensitive 
to any potential LPI/LLI-violating effects. 

The next generation of gravitational experiments will significantly extend our 
current understanding of the empirical foundations of the EEP. A proposed Eotvos 
experiment in space, known as the Satellite Test of the Equivalence Principle (STEP) 
attempts to test WEP to one part in 10^^. The precision of gravitational red shift 
experiments could be improved to one part in 10® by placing a hydrogen maser clock 
on board Solar Probe, a proposed spacecraft (see ref. []^ and references therein). 

The dominant form of energy governing the transitions these experiments probe is 
nuclear electrostatic energy, although violations of WEP/EEP due to other forms of 
energy (virtually all of which are associated with baryonic matter) have also been esti¬ 
mated ip. However there exist many other physical systems, dominated by primarily 
non-baryonic energies, for which the validity of the EEP is comparatively less well 
understood p. Such systems include photons of differing polarization [0, antimat¬ 
ter systems 0 , neutrinos , mesons 0, massive leptons, hypothesized dark matter, 
second and third generation matter, and quantum vaccum energies. Indeed, potential 
violations of the EEP due to vacuum energy shifts, which are peculiarly quantum- 
mechanical in origin [i.e. do not have a classical or semi-classical description) provide 
an interesting empirical regime for gravitation and quantum mechanics. 
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In this paper we investigate the effects that EEP-violating couplings have on 
Lamb-shift transition energies. Such transitions arise solely due to the radiative 
corrections inherent in quantum electrodynamics. A test of the EEP for this form 
of energy therefore provides us with a qualitatively new empirical window of the 
foundations of gravitational theory. 

The Lamb shift is the shift in energy levels of a Hydrogenic atom due to radiative 
corrections. Such energy shifts break the degeneracy between states of with the 
same principal quantum number and total angular momentum, but differing orbital 
and spin angular momenta. The best known example is the energy shift between 
the 2 S'i/ 2 and 2Pi/2 states in a Hydrogen-like atom, which arises due to interactions 
of the electron with the quantum-held-theoretic fluctuations of the electromagnetic 
field. For metric theories, the lowest order contribution for the Lamb shift is 1052 
MHz for hydrogen atoms. There is a 5 MHz discrepancy with the experimental value 
of 1057.845(9) MHz [|^ or 1057.851(2) MHz [|^ , that can be improved with the 
inclusion of higher order terms and corrections coming from the structure and recoil 
of the nucleus. 

Any breakdown of LPI/LLI is determined entirely by the form of the couplings of 
the gravitational held to matter since local, nongravitational test experiments simply 
respond to their external gravitational environment. To explore such effects it is 
necesssary to develop a formalism capable of representing such couplings for as wide a 
class of gravitational theories as possible. We consider in this paper Lagrangian-based 
theories in which the dynamical equations governing the evolution of the gravitational 
and matter helds can be derived from the action principle 


6 J d^x C = S J d^xiCc + Cng) = 0 ■ 


( 1 ) 


The gravitational part Cg of the Lagrangian density contains only gravitational helds; 
it determines the dynamics of the free gravitational held. The nongravitational part 
NG contains both gravitational and matter helds and dehnes the couplings between 
them. The dynamics of matter in an external gravitational held follow from the action 
principle 

6jd^xCNG = 0 ( 2 ) 

by varying all matter helds in an external gravitational environment. 

We work in the context of a wide class of non-metric theories of gravity as de¬ 
scribed by the THefi formalism |^. Phenomenological models of C^g provide a 


general framework for exploring the range of possible couplings of the gravitational 
held to matter and, thus, the range of mechanisms that might conceivably break 


2 





LPI or LLI. The THefi formalism is one such model. It deals with the dynamics of 
charged particles and electromagnetic fields in a static, spherically symmetric gravi¬ 
tational field. In addition to all metric theories of gravitation, the THefi formalism 
encompasses a wide class of non-metric theories. 

A quantum-mechanical extension of the original classical THe^ formalism was 
developed by Will to calculate the energy shifts (due to e.g. hyperhne effects) 
in hydrogenic atoms at rest in a THefi gravitational field. Since the ticking rate of 
a hydrogen-maser clock is governed by the transition between a pair of these atomic 
states, this extension can be used to determine the effect of the gravitational field on 
the ticking rate of such clocks. This provides a basis for a quantitative interpreta¬ 
tion of gravitational redshift experiments which employ hydrogen-maser clocks, for 
example, the gravity probe A rocket-redshift experiment |^. Such experiments are a 
direct test of LPI. 

This formalism was further extended by Gabriel and Haugan who calculated 
the effects the motion of an atomic system through a gravitational field would have 
on the ticking rate of hydrogen-maser and other atomic clocks. Their extension can 
be used to compute energies of hyperhne and other energy shifts of hydrogen atoms 
in motion through a THefi held. Here the physical ehect under consideration is 
time dilation rather than the gravitational redshift. When LLI is broken, the rates 
of clocks of diherent types that move together through the gravitational held are 
slowed by diherent time-dilation factors. This nonuniversal behavior is a characteristic 
symptom of the breakdown of LLI [^, just as nonuniversal gravitational redshift is 


the hallmark of LPI violation 14 


We are concerned in this paper with extending this analysis to the Lamb shift, an 
energy shift whose origin is due to radiative corrections. We compute the Gravitation¬ 
ally Modihed (GM) Lamb Shift in a THefi held, and then discuss experiments which 
could potentially measure such ehects. We hnd both EEP-violating contributions to 
the Lamb shift from the semiclassical THefi Hamiltonian and its radiative correc¬ 
tions. The semiclassical contribution violates LLI only and is isotropic; the radiative 
corrections violate both LLI and LPI and are not isotropic. These contributions are 
functions of non-metric parameters which arise in the leptonic sector of the standard 
model, and so are not constrained by previous high-precision experiments which have 
set stringent bounds for analogous parameters in the baryonic sector [|]. Of course 
all such contributions vanish for metric theories. 

In order to calculate the (GM) radiative corrections, we shall modify the Feynman 
rules of Quantum Electrodynamics (QED) within the context of the THefi formalism. 
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Although we cannot use LPI/LLI symmetries, the gauge invariance of the theory is 
still present. We shall be concerned with the one photon contribution to the (GM) 
Lamb shift up to order ma{Za)‘^, with the nucleus treated as a fixed point charge. 
We do not include further (higher-order) refinements, since we are interested in the 
role of Lamb shift energies in the investigation of possible LPI/LLI violations and so 
expect any such violations to be qualitatively different from higher order corrections. 

Our paper is organized as follows. In Sec. II the THefi action is introduced 
and extended to frames moving with respect to the preferred frame defined by the 
T Hefi gravitational field. This formalism is then used to calculate the electromagnetic 
fields produced by a point-like charge and to formulate (GM)QED. In Sec. Ill the 
(GM) Dirac equation is used to find the energy levels of hydrogenic atoms, and we 
compute the radiative corrections for those states in Sec. IV. In Sec. V the GM 
Lamb shift is related to redshift and time dilation parameters to study possible LPI 
and LLI violations respectively. Final conclusions are presented in Sec. VI. Several 
appendices summarize details of our calculations. 


II (GM) Action 


The THefi formalism was constructed to study electromagnetically interacting charged 
structureless test particles in an external, static, spherically symmetric (SSS) gravi¬ 
tational field, encompassing a wide class of non-metric (and all metric) gravitational 
theories. Originally employed as a computational framework designed to test Schiff’s 
conjecture |^, it permits one to extract quantitative information about the impli¬ 
cations of EEP-violation that can be compared to experiment. It assumes that the 
non-gravitational laws of physics can be derived from an action: 

Sng = - I] / dt (T - HvlY^^ + f dt <dl^«) 

+^Jd^x{eE^-By^i), (3) 


where Ca, and are the rest mass, charge, and world line of particle a, 

= t, = dx^/dt, E = — VAq — dA/dt, B = V x A. The parameters e, and p 

are arbitrary functions of the Newtonian gravitational potential U = GM/r (which 
approaches unity as U —0), as are T and El which in general will depend upon the 
species of particles within the system (leptons in the present case). 

A quantum mechanical extension of the action (|^) which incorporates the Dirac 
Lagrangian was used by Will IM to study the energy levels of hydrogen atoms. In 
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that case a local approximation to the action is employed. The spacetime scale of 
atomic systems allows one to ignore the spatial variations of T, H, e, /i, and evaluate 
them at the center of mass position of the system, X = 0. This work was further 


extended by Gabriel and Haugan who showed that after rescaling coordinates, 
charges, and electromagnetic potentials, the held theoretic extension of the action (^) 
can be written in the form 

1 


S' = y d'^x %lj{i ^ + e 4 ~ + 2 / 




(4) 


where local natural units are used, 4 = and = i^o/^o^oh-o with the subindex 

“0” denoting the functions evaluated at X = 0. The parameter c is the ratio of 
the local speed of light to the limiting speed of the species of massive particle under 
consideration. 

The action (|^) (or (|)) has been widely used in the study of LPI/LLI violating 
effects such as the effect of non-metric gravitational helds on the differential ticking 
rates of different types of atomic clocks, a violation of LPI [Q. An analysis of the elec¬ 
trostatic structure of atoms and nuclei in motion through a THefi gravitational held 
using shows that the non-metric couplings encompassed by the THe^ formalism 

m 


can also break LLI 


This symmetry is broken when the local speed of light 


= (/ioCo)”^^^ dihers from the limiting speed of a given species of massive particle 
Co = (To/hfo )^^^5 the latter being normalized to unity in (^). Further implications 
of the breakdown of LLI on various aspects of atomic and nuclear structure have 
also been investigated. Shifts in energy levels (including the hyperhne splitting) of 
hydrogenic atoms in motion through a THefi gravitational held have been calculated 
pr5| by transforming the representation of the action (^ to a local coordinate system 
in which the atom is initially at rest and then analyzing the atom’s structure in that 
frame. The local coordinate system in which the THefi action is represented by Eq. 
(^), is called the preferred frame; moving frames are those systems of local coordinates 
that move relative to the preferred frame. 

In the present work we generalize this analysis by using the action (^) to study 
radiative corrections to bound state energy levels in hydrogenic atoms. We follow 
the scheme given in Ref. |^|, and analyze the atomic states in moving frames whose 
velocity is u. 

Consider an atom that moves with velocity u relative to the preferred frame. The 
moving frame in which this atom is initially at rest is dehned by means of a standard 
Lorentz transformation. A convenient representation of the THefi action in this 
new coordinate system if the nongravitational helds 'ip, A, E, and B transform via 
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the corresponding Lorentz transformations laws for Dirac, vector, and electromagnetic 
fields is, to O('U^), 


S = j d'^x ^ + e 4 — m)ip + J d‘^x 


+ 

+ 




^ - {u ■ Ef + (1 + u^)B‘^ 


{u-Bf + 2u- {E X B)) 


(5) 


where is the electromagnetic 4-current associated with some external source (taken 
here to be a pointlike spinless nucleus). In our formulation, all non metric effects arise 
from the inequality between cq and c* in the electromagnetic sector of the action. 
The dimensionless parameter ^ = 1 — (c*/co)^ = 1 — measures the degree to which 
LPI/LLI is broken for a given species of particle. The natural scale for ^ in theories 
that break local Lorentz invariance is set by the magnitude of the dimensionless 
Newtonian potential, which empirically is much smaller than unity in places we can 
imagine performing experiments |l|. We are therefore able to compute effects of the 
terms in Eq. (|^) that break local Lorentz invariance via a perturbative analysis about 
the familiar and well-behaved c —^ 1 or ^ 0 limit. 

The fermion sector of the action (^) implies that the equation of motion for the 
'ip held is simply the Dirac equation coupled in the usual fashion to the potential 
Afj,. On the other hand, the pure electromagnetic part of the action is modihed with 
an extra term proportional to the small (species-dependent) parameter This will 
affect the electromagnetic held equations, and the photon propagator. In both cases 
we can calculate ehects of the additional terms perturbatively. 

The held equations coming from the action 


are 
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( 6 ) 


V-E = p + i u-V{u-E)-u-V X B— ifV ■ E , 

W X B — E = j + ^ V X {u X E) + u X V{u ■ B) + {1 + {f)'V x B 
+ u^E — u{u ■ E) — u X B 

where p and j are the charge density and current associated with the fermion held 
plus and external source (such as a nucleus.) Perturbatively solving these equations 
for electromagnetic potentials produced by a pointlike nucleus of charge Ze at rest in 
the moving frame yields 


Ao = [I - + {u ■ hf)](j) = (j) + i(j)' 

A = ^[u + n{u ■ n)\(l) = p^A' 


( 7 ) 
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where n = x/\x\, cf) = Ze/ 47 r|x|, and V ■ A = 0. Note that Eq. (|^ agrees with the 
corresponding resnlt from Ref. 


15 


The primed helds in Eq. (|^) signal a breakdown of LLI. Conseqnently we expect 
that this electromagnetic potential will modify the energy states of hydrogenic atoms 
prior to the inclnsion of radiative corrections. We shall calcnlate these effects for the 
Lamb shift in the next section. In order to hnd the radiative corrections to these 
energy levels we mnst reformnlate Qnantnm Electrodynamics (QED) according to 
the action (^). Since the fermion sector of the action does not change, the fermion 
propagator is nnaltered; only the photon propagator needs to be modihed. 

To hnd the photon propagator, we go back to the action (^) and add a gange 
hxing term of the form 


Sgf — 


(Tx 


{I - i){d ■ Af - Aod ■ A 


after which the resnlting electromagnetic part can be written as 

■1 




EM 


= j d^x ^A^d^d.A^ 

+ ^{A^dod^A>^ + Aod>^d^A^-A,d'^d,A'^) 


( 8 ) 

(9) 


where we have integrated by parts and neglected snrface terms. 

This action is still given in preferred frame coordinates. We can go the moving 
frame by performing the Lorentz transformations 


Ao^ Aq = 7(Ro -u- A) =ji3 ■ a (10) 

do^d'Q = 7(^0 -u-V) = 'yl3 ■ d 


where 7 ^ = 1/(1 — u^) and = (!,«); henceforth = 1 — v?. Transforming Eq. 
d^) by nsing Eq.([T0|) gives 

Sem = II d‘^xA^IC,,A'^ ( 11 ) 

where (in momentnm space) 






( 12 ) 


where 7 ^ 1 , is the Minkowski tensor with a signatnre (H-) and is the inverse 

of the photon propagator G^y. Therefore after solving 




SI 


(13) 
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we find up to first order in ^ 


G,, = -(l + i)^ + i-h 


.7 


/c^ 


'P 


+ Pi^Pu 


(14) 


The terms proportional to ^ in Eq. (^) signal the breakdown of both LPI and 
LLI, since those terms are still present even if u = 0. The (GM) QED then differs 
from standard QED only in the expression for the photon propagator; the fermion 
propagator and Feynman rules are unchanged. 

As the Lamb Shift is the shift between the 2S'i/2 and 2Pi/2 states, and since the 
Dirac equation for a Coulomb potential predicts those states to be degenerate, the 
difference between them in metric theories comes only from radiative corrections. 
For non-metric theories which can be described by the TiLe/i formalism these energy 
levels will be modihed by the EEP-violating terms introduced in the source (Eq. 
®), removing this degeneracy before introducing radiative corrections. Note that 
the fermion sector of the THefi action does not change and therefore neither does 
the Dirac equation. The preferred frame effects appear only in the expression for 
the electromagnetic source produced by the nucleus. We shall now evaluate this 
contribution. 


Ill (GM) Dirac States 

The Dirac equation in the presence of an external electromagnetic held still reads as 
in the metric case: 


H\n) = {a ■ p + {3m — eA° + ea ■ A)\n) = En\n) (15) 

where the various symbols have their usual meaning. 

The (GM) energy levels of hydrogenic atoms are found by solving (|^ in the pres¬ 
ence of the electromagnetic held (J^) produced by the nucleus which entirely accounts 
for the preferred frame ehects. If we replace Eq. (J^) in (^), the Hamiltonian can be 
written as 

H = Ho + H' = + ea-A' (16) 

where Ho corresponds to the standard Hamiltonian (with Goulomb potential only), 
and the primed helds are dehned as in Eq. (^. In terms of the known solutions for 
Ho\ n)^ = E^\n)^, we can perturbatively solve Eq. (|^ by writing 

\n) = |n)° -I- 
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(17) 





with 


K = = K''"'+ £:'"> 

^ir\H'\nf 


W)' = E 

r^n 


E2 - E2 


k) 


(18) 

(19) 


where E'J'^'^ and account for the contributions coming from the respective 

electric and magnetic potentials. 

We now proceed to calculate the energy levels related to the Lamb shift states. To 
obtain these, we find it convenient to use the exact solution for the Dirac spinor In)*^, 
expanding the hnal answer in powers of Za to 0{{Za)^). The relationship between 
this approach and an alternate one in which the Hamiltonian is first expanded in 
powers oi Za using a Foldy-Wouthuysen transformation is discussed in appendix 
The unperturbed Dirac state |?7,)° can be expressed as: 

-iEij{r) a-h \l-,jm) J 

where jl'jj'm) is the spinor harmonic eigenstate of J^,L^ and with respective quan¬ 
tum numbers j, I and m. The functions E and G can be written in terms of confluent 
hypergeometric functions that depend in a non-trivial way on Za for a given I and j 

m- 


Inserting the helds from (^ and (p0|) in E'^, we write 

E'J'^'^ = {Rgg + Rff) {3rn-,l\u^ + {u ■ nf\l-Grn) 

_ —iRQpi^jm]l\{a ■ n){(T ■ u) + u ■ h\l-Gm)+ h..c. 

where “h.c.” means Hermitian conjugate and where 


Rgg = ^ 


G—^Gr^dr 

r 


( 21 ) 

( 22 ) 

(23) 


with Rff and Rgf defined in an analogous manner. 

We now evaluate this energy for the 251/2 and 2Pi/2 states in this semiclassical 
approximation, prior to the inclusion of any radiative corrections. Since the angular 
operator in (^) has odd parity (as given by h), it is straightforward to show that 
the magnetic contribution E'J'^'^ = 0, so E'^ = E'J^^'^ for any state. Using the 
corresponding expressions for the harmonic spinors and the E, G functions in 


for each Lamb state [|^, we hnd 

1 


E' 


2S 


1/2 


7 19 

= -u^m(ZaY l + (-1- ){ZaY 


E' 


6 

1 


'2P1/2 = ^um{Za) 


l + (l5 + i5)(2a.) 


+ 

+ 


(24) 

(25) 
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where we have expanded the exact solutions for Rgg and Rff in powers of (Za)^, 
and kept the hrst relativistic correction only. The angular integration and the Rgg 
term are the same for both states, and so the non-relativistic limit is still degenerate 
for them. However the hrst relativistic correction coming from the Rff factor breaks 
the degeneracy, yielding 

- Bzp.,, = (jm{Zar + O ((Za)») (26) 

We obtain the result that the 2S'i/2-2Pi/2 degeneracy is lifted before radiative 
corrections are introduced. This ‘semiclassical’ nonmetric contribution to the Lamb 
shift is isotropic in the 3-velocity u of the moving frame and vanishes when u = 0. 
Hence it violates LLI but not LPI. 

In order to proceed to a computation of the relevant radiative corrections, we need 
to hnd the perturbative corrections for the energies and spinor states given by ([T^ ) 
and (|^) respectively. The radiative correction 6En to the Dirac energy can be 
formally expressed as 

6En = {n\6H\n) (27) 

where 6H accounts for the loop contributions as given by the gravitationally modihed 
QED. Since EEP violating effects appear in both the photon propagator and the 
classical electromagnetic held, we expect 

6H = 6H^ + ^6H' (28) 


In addition, the state \n) may be analogously expanded. Up to hrst order in we 
can therefore write (^) in the form 


5E^ =° {n\5H^\nf + ^ \^{n\5H'\nf + {\n\5H^\n)' + h.c.} 


(29) 


The contributions from the \n)' states are of the same order of magnitude (in terms 
of powers of Za) as the 5H' terms and so cannot be neglected. This may be seen 
by noting that, apart from the u dependance, 0' ~ 0 and so °(?7,|i7'|r)° ~ - E«. 

Inserting this in (pTl) proves the statement. Note that the ehect of the \n)' states 
was overlooked in Ref. |^. If we identify 5H where P(/i/) represents the 


perturbation to the Dirac Hamiltonian due to the spin of the nucleus, then by the 
same arguments as before, we can show that the term {°(n|i7(/i/)°|n)' -|- h.c.} was 
omitted in the corresponding expression for the hyperhne energy. 
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IV (GM) Radiative Corrections 

To lowest order in QED there are two types of radiative corrections to the energy 
levels of an electron bound in an external electromagnetic potential: the vacuum 
polarization (If) and self-energy (S), along with a counterterm {6C) that subtracts 
the analogous processes for a free electron. These contributions are illustrated in 

Fig.l. 



Figure 1: Radiative corrections of order a : (a) self-energy and (b) vacuum polariza¬ 
tion. 

The energy shift due to these contributions for the state \n) can then be written 
as 

5En = SEs + 6Ep (30) 

where 

6Es = {n\E-6C\n), (31) 

which corresponds to the self-energy contribution in Fig. 1(a) minus the correspond¬ 
ing counterterm, and 

5Ep = (n|n|n), (32) 

which is the vacuum polarization contribution illustrated in Fig. 1(b). 

In Fig. 1 the bold line represents the bound electron propagator. This propagator 
can be written in operator form as (^— Y — m)~^, with 

V^{x) = —eA^{x) and = {En^p) 

where is the external electromagnetic potential. Here E^ is the total energy of 
the state |n), which satishes the Dirac equation (|i— Y — m)\n) = 0 

Fq. ( |30[) represents the one loop correction (one power of a) to the atomic energy 
levels as given by En- We are interested in obtaining the “lowest order” Lamb shift, 
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which is the a{ZaY contribution. (There are still more approximations that come 
after expanding the bound propagator, which introduce additional nonanalytic terms 
in the expression for the Lamb shift that behave like a{Za)‘^ln{Za)). 

The GM radiative corrections are found by evaluating (|30| ) where the external 
electromagnetic potential and the photon propagator are respectively given by Eqs. 
(|^) and (|^. All expressions will be expanded in terms of the LPI/LLI violating 
parameter and the velocity of the moving frame u up to O(^) and as implied 

by© and (p!^ . EEP-violating effects are all contained in the terms proportional to 
these quantities. 

A variety of methods are available for evaluating the corrections in (pOD, each 
differing primarily in the manner in which the bound electron propagator is treated. 
We shall follow the method of Baranger, Bethe and Feynman |]^ (hereafter referred 
to as BBF), in which the corrections in (^) are separated into a term in which the 
external potential acts only once, and another term in which it acts at least twice. 
This latter ‘many-potential’ term can be further separated into a nonrelativistic part, 
and a relativistic part which can be calculated by considering the intermediate states 
as free. This approach is sufficient for the lowest order calculation we consider here. 
We now proceed to outline the main steps of this method. 

The self-energy term in Eq. ( pOD can be written as 

1 


SEc = 


a 


(Ek 


m 


-lAn) - {n\6C\n) 


(33) 


J ^ Y— 

This expression gives a complex result for the level shift, since the denominators in 
the integral each have a small positive imaginary part. The resulting imaginary part 
of 5Es represents the decay rate of the state \n) through photon emission. The Lamb 
shift refers to the real part of the shift, and only that part will be retained in the 
computation of Eq. (|3^) . 

The difficulty in evaluating Eq. (|^ arises entirely from choosing a convenient 
expression for the bound propagator. The integrand in (^) is rearranged in order to 
obtain one part which is of first order in the potential (5Ei), and another part {SE^) 


which contains the potential at least twice. Using the identity |18 

^{,6 + 0 AO + 0 Y 


O = {A- m)- 


pI-pI 


pI-pI 


'a - m), 


(34) 


to re-express 7 ^ and 7 y in (|33|) and respectively identifying pb = p, pa = p — k, and 
Pb = p — k^ Pa = p yields after some manipulation 


5Es = 6 El + 6E2, 


(35) 
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where 


(36) 


= - f d'pd'p'xl;Jp'){h + h + h}i’n(p), 

71 J 

with 


and where 


h 

h 

h 


i 

47|-2 


2p;, y^Py- hp 


r 


kP — 2p'-k — 2p ■ k 

f 2p^ 2p^- 

J k‘^ — 2p ■ k k"^ — 2p ■ k 

- 7m " 


G^''{k)d^k 

G^'^{k)<fk 


dTT^ J k? — 2p ■ k 


j^Gf^’'{k)d^k - SG 


iSSj = ^/ i-^{f)M^(p',p' -s'-k) 

X K^{Eo — ko-,p'— s'— k,p + s —k) 

X MI{p + s — k,p)'ipn{i^G^’^{k)d'^kd^pd^p'd^sd^s' 

= {mkXm) 


(37) 


(38) 


with 


Mf,{p',p-k) 

Ml{p' - k,p) 


np'-p) 

np'-p) 


2Pm - 7m 

2p ■ k — k? 

2py- ji'yy 

k"^ — 2p ■ k 


2p), - 7 m 

2p' ■ k — k^ 

‘^P'u- ^Iv 

k? — 2p' ■ k 


f{p'-p) 

f{p'-p) 


The quantity is dehned as —iK^ = y — m) where in momentum space 

K'l = 5(E' - E)KX(E-,p\^. 

In Eqs.(PB[) and ( P8|) the p’s have time component En and the s’s have time 
component 0. Note that the above derivations are independent of the specihc form 
of the photon propagator 

Further evaluation entails a lengthy computation which in principle is analogous 
to that of BBF. In practice though, the calculation is substantially more complicated 
than in the metric case due to the additional non metric terms present in the pho¬ 
ton propagator and the electromagnetic source related to a charged point particle. 
Regularization and renormalization procedures have to be modihed accordingly as 
well. Details involving the subsequent computation of the self energy (and vacuum 
polarization) term are given in appendix 
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The final result for the loop corrections related to the Lamb shift is of the form 


= 6E. 


2S 


1/2 


-6E. 


2Pi 


1/2 


(39) 


where each term is obtained from eq. 


(and its relevant subsidiary equations) as 


calculated for the corresponding atomic state. By adding the “semiclasical” correction 
coming from the Dirac level (labeled by {D) in Sec. Ill), the total Lamb shift reads 


+ ^EiQ) 


AEr, = AE 


2.084 + In — + e 




4.534 + ? In At 




(40) 


+ r 


- - 3.486 + ? In At - 0.011 cos^ 9 


a 




+ UiUjAe^^ I 


where we have introduced the dimensionless parameter Ae*-^ = 2AE^^/{{Za)^m^) (see 
), and used Eqs. (|B47]) and ([B48|) in the evaluation of (^) through Eq. 


The former result is the energy shift associated with the particular states in 
However in Eq. (|B43| ) we have derived a general expression for the one-loop radiative 
corrections related to any atomic state. These are 


TZ. 4 {Za^a 
dEno = - T — m 


dvr 




m 


—- —+ (l + -01n(- 
30 30 ^ 2^^ ^2EI‘0' 


+ 0{u^ 


for I = 0, and 


4 {Zafa 

dEni = - ^—m 


Ott 


n'^ 


3^,, ,Z‘^Ryd^ ^ 3 Cij ^ ^ 


(1 + 2^)ln( ) + 




(41) 


(42) 


for / 7 ^ 0; where we have not explicilty written the terms proportional to the moving 
frame velocity. Here 


Cu = 


l/(/-|-l) for j =/-I-1/2 
— 1// for j = I — 1/2 


(43) 


and E* is defined by 


Values for this reference energy can be obtained from 
Ref. [T^ up to states with n = 4. 

Note that in addition to the explicit dependence on the frame velocity in Eq. (|40D , 
there exists a position dependence hidden by the rescaling of the original action (Eq. 
d^)), which was considered locally constant throughout the computation. The full 
THey parameter dependence in Eq. (|40|) can be recovered by replacing 


a 


"eV T’ 


m 


m 




AEt 




(44) 
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in the preceding equations. 

Note that ^ in Eq. (^) accounts for any EEP violation coming from a non- 
universal gravitational coupling between photons and leptons. A further distinction 
can still be made between leptons and antileptons. In principle a matter/antimatter 
violation of the EEP could be measured in a Lamb shift transition, through the 
appearance of virtual positron/electron pairs in the vacuum polarization loop contri¬ 
bution 1^. This will add a non metric term to Eq. (^OD, of the form (see appendix 
y for more details): 


AB<+' = 


-c 


e+ 


m 

UOn 


{ZaYa{l + 2|-up) 


(45) 


where = 1 ~ Ce_/ce^ accounts for the difference between the limiting speed of 
electrons (ce_ = Cq) and positrons (ce_,_). 

We turn next to the question of relating the Lamb shift to observable quantities 
in order to parameterize possible violations, of the EEP. 


V Test for LPI/LLI Violations 


We begin by considering a general idealized composite body made up of structureless 
test particles that interact by some nongravitational force to form a bound system. 
The conserved energy function of the body E is assumed to have the quasi-Newtonian 


form 16 


E = MkcI - MxU(X) + -Mx\Vf + .. 


(46) 


where X and V are respectively the quasi-Newtonian coordinates and velocity of the 
center of mass of the body, Mr is the rest energy of the body and U is the external 
gravitational potential. Potential violations of the EEP arise when the rest energy 
Mr has the form 

MrcI = MocI-Eb{X,V) (47) 


where Mq is the sum of the rest masses of the structureless constituent particles and 
Er is the binding energy of the body. It is the position and velocity dependence of 
Er, which signals the breakdown of the EEP. Expanding Er in powers of U and 
to an order consistent with (|4^) we have 

Er{X, V)=El + (48) 

where is the external gravitational potential tensor, satisfying 7/** = U. The 
quantities Snip and Srri/ are respectively called the anomalous passive gravitational 


15 







and inertial mass tensors. They depend npon the detailed internal structnre of the 
composite body. In an atomic system they can be expected to consist of terms 
proportional to the electrostatic, hyperhne, Lamb shift, and other contribntions to 
the binding energy of an atomic state. 

In a gravitational redshift experiment one compares the local energies at emission 
Ef>m and at reception Erec of a photon transmitted between observers at different 
points in an external gravitational held. The measured redshift is dehned as 

Erec 

Eern 


Using (^61) (with U = 0) to relate the transition energies at the two different points, 
this parameter can be expressed as 


Z = Af/(l-S), 


5m% Af/o 
A^o AU 


(49) 


Clearly Z depends (through Srnpj upon the specihc test system used in the experi¬ 
ment. An absence of LPI violations will mean H = 0, and so Z will be independent 
of the detailed physics underlying the energy transition . 

The LLI violations may be empirically probed through time dilation experiments. 
These experiments compare atomic energy transitions as measured by the moving 
frame {AEb) and preferred frame {AE^), which can be related via 

A£;b = AB»(i-|>1-1]U) (50) 

with the time dilation coefficient A dehned by 

^ 6mf 

^ ^ “ AE% 1/2 

Here represents the diherence between the anomalous inertial tensors re¬ 
lated to the atomic states involved in the transition. The coefficient A represents 
the dilation of the rate of a moving atomic clock whose frequency is governed by the 
transition. Since the anomalous mass tensor is not isotropic, A depends upon the ori¬ 
entation of the atom’s quantization axis relative to its velocity through the preferred 
frame. Note that if LLI is valid the anomalous inertial mass tensor associated with 
every atomic state vanishes, so that A = 1. 

Here we consider the possibility of employing the Lamb shift as the atomic tran¬ 
sition governing the appropriate experiment. To do so we must compute the relevant 
S and A coefficients respectively. 
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In order to calculate the corresponding related to the Lamb shift, we must 
hnd the manner in which AE^ varies as the location of the atom is changed. Setting 


n = 0 in (|40D and performing the rescaling given in (^4]), we obtain 

T 


A-Bi = ClA |i + + 4(1 + 

with Sl = ^{Za)^a/b, and 


In e 


H 


(52) 


a = 6(—4.534 + - In b = l/(—2.084 + In 

2 a‘‘ 

where Sl represents the metric value (within the given approximations) for the Lamb 
shift. Note that there is still a position dependence in (|52[) through the dehnition of 


I lie 


(53) 


We recall that the total energy of the system can be expressed in term of 

E = mVf + AEl + ■■■ (54) 

where ■ ■ ■ represents other contributions for the binding energy of the system. 

The functions T, H, e and /r, considered to be functions of U and evaluated at 
the instantaneous center of mass location X = 0 for purposes of the calculation of 
AEl, are now expanded in the form 


T(U) = r„ + r'g„. X + 0(m ■ xf 


(55) 


where go = Vt/|jY=o’ = ^lx=o) — dT/dU\ji^Q. It is useful to redehne the 

gravitational potential U by 

1 Ti 

(56) 


1 T' 


whose gradient yields the test-body acceleration g. 
If the above is used to expand (Q), we get 


E = {m + £^l)(1 -U) + - a- 26)ro - aAoj 


(57) 


where we have used (^4|); and neglected terms proportional to since the main 
position dependence parameterization is given in terms of: 


° u ^eo 2To 2770^’ 


A = ^^0 / ho I ^(3 _ H'o 1 

U^li, 2To 2Ho' 


(58) 
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If we now identify (^Tf ) with Eqs. (|i6|) and (|i8|), we can obtain the corresponding 
Lamb shift contribntions to the binding energy and anomalons passive mass tensor 
as 


= -£l (59) 

— a — 26 )ro — aAoj 


This resnlt was hrst presented in Ref. , where in (^) we have corrected the latter 
for a sign error in the coefficient mnltiplying Aq and a missing factor b in the Tq term. 
Inserting (p9D in (|i9|), we obtain 


= 3.424 To - 1.318 Ao 


( 60 ) 


as the LPI violating parameter associated with the Lamb shift transition. Note that 
if LPI is valid then Pq = Aq = 0. 

In comparing the resnlt (|60D to anomalons redshift parameters compnted for other 
systems, it is important to note that we are working with nnits that are species 
dependent. Recall that the choice of cq = 1, and the redehnition of the gravitational 
potential (|5B[ ) involves the T and H fnnctions associated with electrons (or more 
generally a given species of lepton). 

Consider, for example, hyperhne transitions (maser clocks). In this case the lep- 
tonic and baryonic gravitational parameters appear simnltaneonsly. This atomic split¬ 
ting comes from the interaction between the magnetic moments of the electron and 
proton (nnclens). The proton metric appears only in the latter, and so it does not 
affect the principal and hne structure atomic energy levels. It is simple to check that 
the hyperhne splitting scales as 


^Ehf = 


£hf 


Hb To ef] 


(61) 


where the label B is added to distinguish baryonic related functions from leptonic 
ones; and £hf depends only on atomic parameters. 

In expanding (^) according to (|5^) , we obtain 


= 4/(1 - Ub) + SkiUB’S.'^’ 


( 62 ) 


with 

= 3Tb-Ab + A 


(63) 
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where Ub, and Kb are the baryonic analogues of (^), and (|5^) respectively. In 
(^) we rescaled the atomic parameters to absorb the THefi functions and chose units 
such that cb = 1. The quantity A is given by 


A = 2 


Za 

T' 

-‘-B 


2(^ 

^Hb 


Ho' 


Ik + Io 

Tb To 


(64) 


and would vanish under the assumption that the leptonic and baryonic THefi parameters 
were the same. 

Turning next to experiments which test LLI, we need to obtain the tensor 6rnj 
appropriate to the Lamb shift. This tensor is obtained after taking partial derivatives 
of AEl with respect to Ui and Uj (note V = u). Substituting the result into (|^ ) 
yields 

iL i 


1-A^ = 


7.757 


- + 3.074 - 0.011 cos^ 0 + ^Ae*^' 

a i/2 


(65) 


for the Lamb shift time dilation coefficient, where Q is the angle between the atom’s 
quantization axis and its velocity with respect to the preferred frame. 

Note that the coefficient Al depends upon Ae®-^, the evaluation of which involves 
the computation of an inhnite sum as given by (|B44|) . The dominant contribution in 


Eq. ([65| ) comes from the Dirac part of the energy (proportional to ^ ), which produces 
an overall shift only. Non-isotropic effects arise solely due to radiative corrections. 

In general, an experimental test of LLI involves a search for the effects of motion 
relative to a preferred frame such as the rest frame of the cosmic microwave back¬ 
ground. A detailed analysis about the interpretation of LLI violating experiments is 
presented in Ref. II. which analyzed experiments concerned with hyperhne transi¬ 
tions, obtaining an expression for the time dilation parameter corresponding to that 
kind of transition^. This parameter is negligible in comparison with other sources of 
energy, such as nuclear electrostatic energy in the case of the ^Be"*" clock experiment 


In summary, we have been able to parameterize EEP violations arising from Lamb 
shift transitions associated with redshift and time dilation experiments. In these types 
of EEP violating experiments one typically looks for variations of the energy shift 
due to changes in either the gravitational potential or the direction of the preferred 
frame velocity. The feasibility of such experiments is hindered by the present level of 
precision of Lamb shift transitions (one part in 10®) in comparison to the magnitudes 
of such changes. In the hrst case, any Earth based experiments will be limited by 

^Note that the expression given there for is incomplete according to discussion presented in 
Sec. Ill 
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the small size of the Earth’s gravitational potential (~ 10“®), which is well beyond 
any foreseeable improvement in Lamb shift precision. Similar problems appear in 
the second case, where the known upper bound |m| < 10 “^ |]l| for the preferred frame 
velocity, leaves no room for any improvement on the EEP violating parameter since 
anisotropic effects go as Clnp. 

However useful information can still be extracted from Eq. (^OD if we use the 
current level of discrepancy between the experimental result and the theoretical 
(metric) value ||^ to bound the nonmetric contributions for the Lamb shift. This 
constrains ^ < 1(1) x 10“®. Similar bounds can be obtained by considering empirical 
information about other atomic states. In this context, the indirect measurement 
of the 15 Lamb shift gives a limit ^ < 1.4(1) x 10“®, and the measurement 
of the 25 i/2 — 2 P 3/2 fine structure interval [^]: ^ < 0.7(1.4) x 10“®. If we drop the 


assumption that positrons and electrons have equivalent couplings to the gravitational 
field [^, we find that there is an additional contribution to (|40|) due to ^e+ 7^ Ce-- 
This contribution arises entirely from radiative corrections and is given by eq. ([451) . 
Making the same comparisons as above, we find the most stringent bound on this 
quantity to be |^e+| < 10“^. 

The previous bounds were obtained by using (|21|) and (|4I|) or (|42|) to calculate the 
corresponding nonmetric Dirac and radiative corrections contributions respectively. 
The 15 Lamb shift experiment, actually measures the transition: {E^s — P 25 ) ~ 
1(^25 ~ -^ 15 )) cind so we use this one to make the comparison, where experimental 
and theoretical values are given in ref.[^]. In the other experiment we need to use 
the non metric part of P 251/2 ~ (= ^^)) namely: 


= ^{ZaYm 


u 


±—(- cos^ 6 

60^2 


1 ) + 0 ((Z«) V) + ^^( 10.434 + 0 (n®)) 

DTT 


( 66 ) 


where the hrst term comes from the Dirac contributions (here + and - label the 
transition coming from the 2 P 3/2 state with \M\ = 3/2 and \M\ = 1/2 respectively) 
and the second one from radiative corrections. Note that the leading anisotropic 
effects stem from the non relativistic contributions, and so their ratio with the metric 
value, 0{m{ZaY), is j{Za)"^) , instead of O(^u^) as for the classical Lamb shift. 
Time dilation experiments will look for changes on the P 


2S 


1/2 


-P. 


2Pc 


3/2 


splitting as the 

Earth rotates, which would single out only the preferred frame contributions. Current 
experiments |23] measure a value of 9911.200(12) MHz for that transition, which gives 


a nominal bound (coming from the experimental error) of cos^ 0 < 1 x 10 ^ for 
the preferred frame part. This bound should improve once appropriate experiments 
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are carried out, since these will look for periodic behaviour which can be isolated and 
measured with high precision. 

Note that an empirical value for the Lamb shift is obtained from Ref.||2^ by sub¬ 
tracting the theoretical result of the hne splitting 2Pi/2 ~ ‘^Psi 2 - Now by following the 
previous formalism we can parameterize the LPI violation in the former experimental 
result through: 




2‘S'i/2“2P3/2 


= {£, + £l)(1 -U)+ U{£,EI + 


(67) 


where we have added the corresponding parameters related to the fine transition 
[|I|: £f and Constraining the ratio of this quantity to a direct measurement of 
the Lamb shift to lie within experimental/theoretical error, we obtain the bound 
—S-f)! = |17(0.576ro-|-1.318Ao)| < 10“®. This result is sensitive to the absolute 
value of the total local gravitational potential whose magnitude has recently 

been estimated to be as large as 3 x 10“® due to the local supercluster ig. Hence 
measurements of this type can provide us with empirical information sensitive to 
radiative corrections that constrains the allowed regions of (Lq, Aq) parameter space. 
Unfortunately the present level of precision in measuring the Lamb shift allows only 
a rather weak constraint. 


VI Discussion 

We have computed for the hrst time radiative corrections to a physical process, namely 
the energy shift between two hydrogenic energy levels that are semi-classically degen¬ 
erate, within the context of the THefi formalism. The corresponding (GM) QED 
was derived, and the (GM) expressions for the propagators were obtained. The non¬ 
metric aspects of a theory describable by the TiLe/i formalism can be all included 
in the photon propagator, given an appropriate choice of coordinates, leaving the 
fermion propagator unchanged. The addition of more parameters to the theory (by 
the THefi functions) entail new renormalizations, where not only charge and mass 
need to be redefined but also the THe^ parameters. 

The approach we took to solve for the semi-classical Dirac energies (sec. Ill) differs 
from the one given in Ref. [|^, in which the Dirac Hamiltonian was expanded us¬ 
ing Foldy-Wouthuysen transformations yielding the first relativistic correction to the 
Schrodinger Hamiltonian (as introduced for example, for the Darwin and spin-orbit 
terms), and subsequently the energies. Instead we began from the fully relativistic 
expression, where the perturbations come only from the preferred frame terms of the 
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electromagnetic potential. Our approach involved evaluating expectation values with 
respect to the relativistic spinors instead of their nonrelativistic extensions (or Pauli 
states). The effects of relativistic corrections such as spin-orbit coupling are therefore 
included exactly in this approach. Once this is done, the hnal result is expanded to 
keep it within the desired order. The semi-relativistic approach is not suitable when 
preferred frame effects are studied. 

Qualitatively new information on the validity of the EEP will be obtained by 
setting new empirical bounds on the parameters and which are associated 

with purely leptonic matter. Relatively little is known about empirical limits on 
EEP-violation in this sector 0. Previous experiments have set the limits [Q = 
|1 — Cg| < 6 X 10“^^ where cb is the ratio of the limiting speed of baryonic matter to 
the speed of light. In our case we obtain an analogous bound on ^ for electrons from 
the difference between current experimental and theoretical values, giving |Q < 10“®. 
Although much weaker than the bounds on ^b, it is comparable to that noted in 


a different context by Greene et. al. 


They considered a similar formalism 


{THefi with u = 0) for analyzing the measurement of the photon wavelength emitted 
in a transition where a mass Am is converted into electromagnetic radiation, thereby 
providing an empirical relationship between the limiting speed of massive particles 
(electrons) and light. 

The breakdown of LPI for the Lamb shift in the context of a nonmetric theory 
of gravity describable by the THefi formalism is embodied in the the anomalous 
gravitational redshift parameter (|60D. Recall that H depends on the nature of the 
atomic transition through the evaluation of the anomalous passive tensor. This tensor 
will have differing expressions for differing types of atomic transitions [Q. An atomic 
clock based on the Lamb shift transition will, in a non-metric theory, exhibit a ticking 
rate that is dependent upon the location of the spacetime frame of reference and that 
differs from frequencies of clocks of differing composition. For example, the gravity 
probe A experiment employed hydrogen-maser clocks, and was able to constrain 
the corresponding LPI violating parameter related to hyperhne transitions: 


This experiment involves interactions between nuclei and electrons and so does not 
(at least to the leading order to which we work) probe the leptonic sector in the 
manner that Lamb-shift experiments would. In general Eq. (^9]) will describe the 
gravitational redshift of a photon emitted due to a given transition in a hydrogenic 
atom; for a hyperhne transition the redshift parameter is (|68D, whereas it is (|60|) for 
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the Lamb shift transition. 

An analogous experiment to test for LPI violations based on Lamb shift tran¬ 
sition energies poses a formidable experimental challenge because of the intrinsic 
uncertainties of excited states of Hydrogenic atoms. Setting empirical bounds on 
by precisely comparing two identical Lamb shift transitions at different points in a 
gravitational potential would appear unfeasible since the anticipated redshift in the 
background potential of the earth 10“®) is much smaller than any foreseeable im¬ 
provement in the precision of Lamb-shift transition measurements ||^. One would 
at least need to perform the experiment in a stronger gravitational field (such as on a 
satellite in close solar orbit) with 1-2 orders-of-magnitude improvement in precision. 
A ‘clock-comparison’ type of experiment between a ‘Lamb-shift clock’ and some other 
atomic frequency standard [|^ is, in principle, sensitive to the absolute value of the 
total local gravitational potential [^, as noted earlier. With this interpretation, 
comparitive transition measurements of the type discussed in the previous section can 
more effectively constrain the allowed regions of (ro,Ao) parameter space than can 
measurements which depend upon changes in the gravitational potential. Of course 
exploiting anticipated improvements in precision of measurements of atomic vacuum 
energy shifts 0 will yield better bounds on ^g- and .^e+ via (^Ol). 

Violations of LLI single out a preferred frame of reference. In fact, the search 
for a preferred direction motivated the most precise tests of LLI performed so far 
P, We have extended the analysis of the effects of motion relative to a preferred 
frame to account for the radiative correction for the atomic energies associated with 
the Lamb shift, as embodied in the expression ([651) . This non-universality reflects 
the breakdown of spatial isotropy for quantum-mechanical vacuum energies. The 
coefficient depends upon Ae*-^, the evaluation of which involves the numerical 
computation of the sum in (|B44|) . Unfortunately, the intrinsic linewidths of the 
relevant states render direct measurement of such effects unfeasible. More precise 
empirical information on the value of ^ can be obtained by precisely measuring changes 
in the T^ 2 Si /2 ~ splitting as functions of terrestrial or solar motions. However 

these effects are insensitive to radiative corrections, depending instead upon the semi- 
classical non-metric effects discussed in section III. 

Finally, we note that our formalism could be applied for muonic atoms. For a 
muon-proton bound system, we will obtain an expression similar to that of ([B43|) , 
but where all parameters refer to muons. For an anti-muon electron bound system 
(a muonic atom) a similar analysis would apply. However in both cases the mass and 
spin of the muon could not be neglected. 
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We expect that the intrinsically quantum-mechanical character of the radiative 
corrections will motivate the development of new LPI/LLI experiments based on the 
Lamb shift transition. In so doing we will extend our understanding of the validity 
of the equivalence principle into the regime of quantum-held theory. 
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Appendices 

A Semi-Relativistic Calculation of Hydrogenic En¬ 
ergy Levels 

Consider a hydrogenic atom immersed in an external gravitational held, moving with 
velocity u relative to the preferred frame. In Sec. Ill we follow a fully relativistic 
approach to solve for the atomic energy levels. That is, we perturbatively solve the 
Dirac equation in the presence of the electromagnetic held of the nucleus, where the 
unperturbed states correspond to the Dirac solution in the presence of a Coulomb 
potential only (the metric case). 

We consider here the use of the Foldy-Wouthuysen transformation in solving (pISl). 
In this approach, we write 

H = Hc + Hmag + H^v + Hso + H£, (Al) 

with 

—O 

P 

He = m + - - eAo 

2m 

' A + A ■ ■ B 
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4m^ 
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Exp 


(A2) 


where is given by Eq. (0)- 

As shown in section III, we can take Hmag 0, since the magnetic held does 
not contribute to the atomic energy levels. We can then group the terms in the 
Hamiltonian as 


H = H^ + Hf 

Hf = Hmv + Hso + Ho (A3) 

where we have dehned the hne contribution to the Hamiltonian (hi/), in order to 
account for the hrst relativistic correction 0{{ZaY) to the atomic energy levels. 

We start writting a formal solution for H\n) = En\n), in term of its non-relativistic 
limit: 

H,\n), = E^Jn),, (A4) 


as 


\n) = \n)c+\n)f, E^ = E^ + c{n\Hf\n)c (A5) 

where the index “/” accounts for the hrst relativistic correction to the states and 
energies. 

Since Aq = 0 + ^0', and so He = H'^ + we do not know the exact solution for 
(|A4:|) , but the perturbative expansion: 

\n)c = |n)c° + i\n)'e E^^ = E<^^ +° {n\H[\n)l (A6) 


where 


P 


Hc\n)e = ^ “ e0)|n)° = 


(A7) 

If we use ([A^) along with Hf = + in (|A5|) , we can hnally write up to O(^), 

E„ = + (A8) 

+ e [“{n| (HJ + H/) |n)» + { "MHpn)', + h.c.}] + 0{(Zaf ) 


We see then that under this semi-relativistic approach, we must address the prob¬ 
lem of hnding the states |n)(,, whose contribution to (A8|) is between the brace brack¬ 
ets. This is equivalent to include the hrst relativistic correction coming after solving 

= (if° + ifj + ---)|n)° (A9) 
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as 


since, we can show 

{c{n\Hf\ny^ + h-c.} = { + h.c.} 

This relation allows us to rewrite part of (|A 8 |) as 


(AlO) 


(All) 


En — (c(^l + /(^l + ■ ■ ■) (-f^c + -^/ + ■ ■ ■) (l’^)c + 1^)/ + ■ ■ ■) 
= ^{n\H'\n)^. 


(A 12 ) 


It is clear then that if we start with the exact solution for the Dirac equation in 
the presence of a Coulomb potential, we can avoid working with the states \n)y Note 
that since we are interested only in the hrst relativistic correction, the result ( |A12|) 
must be expanded to 0{{Za)‘^). 

Unfortunately for hyperhne or Lamb shift energies, the effect of the primed states 
cannot be removed, since they both come from perturbations to the (known) rela¬ 
tivistic solution of the Dirac equation in the presence of a Coulomb potential only. 

A semi-relativistic expression for the Hamiltonian of a hydrogenic system was 


worked out in Ref. |T^, where the effects of nuclear spin (hyperhne effect) were also 
included within the context of LLI violations. The result presented there for the 
atomic energy levels is incomplete though, since the contribution of the prime states 
was overlooked, as discussed at the end of Sec. III. 


B Loop calculations 

Given the form of the photon propagator (JT^), it is convenient to divide the calculation 
into two parts 

6Es = 6Ef^ + 6Ef'^ (Bl) 

where 5E^^'^ groups the contributions of the terms proportional to in whereas 
SE^^^ contains those proportional to 7 ^ = 1/(1 — u'^) and We are interested in 
solving for the shift in energy levels up to hrst order in so it is enough to consider a 
Coulomb potential as the source for part B, while for part A the full source as dehned 
in Eq. (|^) needs to be included. 

We mention again that we are interested in calculating the CM Lamb shift to 
lowest nontrivial order in a, i.e. up to 0{a{Za)^). To this order, we can use the 
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nonrelativistic expressions for both the large and small component of the electron 
spinor So for example, if we make the snbstitntion 


= {Zam) 


(B2) 


where is a dimensionless spinor whose hrst two components are of order nnity, 


and the last two are of order Za, we can assign orders to the various terms according 
to 


Pi ~ Zam, Eq — m {Za)‘^m 
cAq d^p' ~ eAi d^p' ~ {Za^m 
-0 li i’nd^p rs-/ Zam. 


(B3) 


These approximations will be used in the sequel to simplify the expressions we obtain. 

B.I Type A Contributions to the Self-energy 

Here we will consider 



(B4) 


and with A^ given by Eq.(|^). This part of the calculation is almost 

identical to that of BBF[^; the only difference is that now we have to consider a 
source that contains a magnetic part in addition to the electric one. 

We begin by computing 6 Ei. Relating the counterterm 6 C to the renormalization 
of the electron mass and regularizing the photon propagator via 



(B5) 


we hnd that I 2 and I 3 in (|36|) become 


I 2 =(1 + 0 -ln(AV/)} 

h =(1 + 0 ^{f(ln(AVpO + ^)+"iln("iVp^)} 


(B6) 


On the other hand, we obtain for Ii 


^ Q 1 /*! rJ'v 1 /*1 

h = “ 3 ^ “ 2 ^'' ^ ^ /o ^ + 4 ^ H^VpI) 



(B7) 


2V -pil-x) -2V ■ px + V ■ p^ 
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where px = xp' + (1 — x)p. 

We can simplify this expression by letting the momentnm operators and 
respectively act on the spinors V’(p') and using the Dirac equation and ([B3[) to 
keep terms up to the desired order. 

Adding together Ji, J 2 , and J 3 we obtain a result correct to order 


iisp = “(1 + 0/v>(p') 


a 


TT 


r- 




(1 + 0(’^l 


-3V^ + 5V^ 

4m 


1 m 1 

3‘"<7’”8 


n), 




i/jif^ap mp 


(B 8 ) 


with q = p' — p, and Note that the term proportional to in Eq. 

(|B 8 |) needs to be evaluated with only the large component of -0 and ~ Vq (70 ~ 1). 

We point out that the initial ultraviolet divergence in (|B^ ) is cancelled after the 
addition of the J’s in (p 8 |). The remaining infrared divergence will be cancelled by a 
similar term which comes from the many-potential part of the level shift. A similar 
cancellation occurs in the non-gauge invariant term present in Eq. (|B^) . These 
cancellations are non-trivial, and provide useful cross checks to our calculation. 

Consider next the evaluation of 6 E 2 . Since the operator satishes the transver- 
sality condition 

k-M = k-M^ = 0 (B9) 

we can write Mq = k ■ M/k^. 

Using 

r = 2U ■ fc 7 ^ - 2V^ j^+ ^ 7 m f, (BIO) 

in the hrst term of Eq.(^) the operator M^can be decomposed into 


with 




(Bll) 

mb = 

r 2p'^ 

k^ — 2p ■ k k"^ — 2p' ■ k 

(B12) 

+ 

^7/.(a,2 _ 2p/ . A; k‘^-2p-k^^^ 


= 

2{y^ f^-V-k^^)/{e-2p-k), 

(B13) 


MB . A; = ■ A; = 0. 


(B14) 


each of which still satishes 












In terms of these operators we now have 
SE 2 = (B15) 


where each term represents a contribution to Eq. (p 8 D involving the products of 
only or or cross terms operators. The simplihcation of these terms is quite 
analogous to that shown in BBF [|^ . The decomposition of the M operator in ( [Bll| ) 
allows one to use simpler expressions for the bound propagator K^. In appendix |B. V| 
it is shown that only in the part {M^K^M^) will it be necessary to use the bound 
electron propagator; in all other contributions it is sufficient to replace by the 
propagator for free electrons, K^. Moreover the main contribution to {M^K^M^) 
arises from intermediate states of the electron with nonrelativistic energy so that both 
and can be replaced by their simpler nonrelativistic approximations. It is 
also shown that the cross term in Eq. ( |B15|) gives a contribution of order a{Za)^ and 
is therefore not relevant in our calculation. According to the above considerations we 
can then approximate Eq. (|B15|) by 


(B16) 


SE^ ~ + {M"KlM") = (M‘) + (M"). 


We start evaluating the hrst term of Eq. ( |B16D . The nonrelativistic prescription 
for KX is given by 


KIr{,x',x) = 
or in momentum space 


Er exp{-iEr{t' - t)) 

0 


for (T — f) > 0 
for it' — f) < 0 


(B17) 


=-i^iPr{p')vl{:^{Er - En +ko) ^ (B18) 

r 

where pr represents the large component of the Dirac spinor. 

In the same nonrelativistic approach reduces to 

m;<™I ~ (p; - p„)ihh^ = (bi9) 

where we have approximated ^17 ~ Vq, because although the magnetic and electric 
potential have the same order of magnitude (as powers of 2 ’a), the 7 matrix mixes 
large components of the intermediate states with small ones and therefore introduces 
corrections one order higher in Za. 
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Therefore, after replacing Eq. 


and 


in Eq. 


we obtain 


(M') = fd‘kG'“'{k) X x: (’yxl’-)(r|fi>> (B20) 

47T % J j, ^0 

where we have neglected the contribution of the photon momentum k to the momen¬ 
tum of the intermediate electron states. This is equivalent to leaving out the factor 
exp(i/c-a;) in the spatial integration. This can be done because k ~ En—Er ~ m^Za)"^, 
which is small compared with the electron momentum p ~ mZa for nonrelativistic 
states. 

Inserting ( |B4|) into ( B20 ), and using Eq. ( B14 ) to relate the temporal component 
of R with its spatial components, which satisfy 


(n|.R|r) = —^{En - E^)(n|^r), 
772/Cq 


we hnd, after integration 


2 a 


{M') = ^(1 + OE - -B«) 


ln( 


/i 


2\En. — Er 


r) + 


6 


(B21) 


(B22) 


where all the states and energies represent the non relativistic limit of the Dirac 
solution. 

Eq. ([B22|) can be simplified by using 


E - B„) = -{n|W„|n), 


(B23) 


which hnally gives 


(M^) = 


a 




(1 + 0 


(ln(^) + b(n|VV„|,i) + C 


with 


(7 = C^\ = 2^(r|p»{n|p,|r)(E, - E^) In 




En — Er 


(B24) 


(B25) 


where E* is a reference energy to be dehned, and has been introduced for later 
convenience. To obtain this result we have neglected the imaginary part of (M^) 
retaining only the leading terms of {M^) in the limit /r ^ 0. 

In computing we can take to be the free electron propagator, which is 


K]^{En - ko]p' - s' - k,p + s- k) = 


— p' + p + s) 
f— ^ — m 


(B26) 
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where 

r^ = (m, s*), s^=p'-s'=p + s (B27) 

upon which (M^^) becomes 

(M^^) = ^ y d^p'd^pd?s^'il){f)Va{f - s^)Np{p^, s^)V^{s^ - (B28) 

with 


N'^{P* 


.) = -- 
i 


^ ^ + m){p0 ^ 

{k'^ — 2p* • ky{k‘^ — 2r-k — sl) 


Gf,^{k)d^k. (B29) 


In the nonrelativistic domain / d^pVa ~ {ZaYm and so the constant value of 
(independent of the momentum and energy of the intermediate states) will already 
yield an overall contribution to Eq. (P28|) of the desired order a{ZaY. Note that 


can be expanded in powers of the momentum p', p or s*, which are of order mZa, 
and therefore any contribution beyond the constant, Za-independent term will be 
of higher order. The same argument can be used to neglect the binding energy of 
the intermediate states. We can therefore evaluate (|B29| ) by approximating p ~ p* 
and p' ~ p* in the denominator of and respectively, so that p* ~ (m, 0) and 
s* 0. 


Evaluating N as in reference [^| we hnd that ( B28 ) becomes 


=-(1 + 


(B30) 


Note that this term will exactly cancel the non-gauge invariant term present in Eq. 


Finally we add Eq. (|B24D to Eq.( |B30|) to obtain 5E^\ and then add it to Eq. 
to give the hnal result for the type-A contribution to the self-energy: 


= 


a 


(1 + f) [c + (In(^) + ^) {n| vv„|n) 


11s 


24 ^ 


37im? 

5m j iip')ia^V^g‘'i’„{p)(fp'd^p 


(B31) 


Apart from the constant (1 factor, there is no formal difference between the 
result ([B31|) for this contribution to the level shift and the standard one . However 


there are implicit differences which appear in the expression for and the solution 
for the Dirac states \n) (in the non-relativistic approach here) in the presence of that 
source. 
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B.II Type B Contributions to the Self-energy 

To solve the type-B contributions we have to consider the photon progator 


= e 


11 




(B32) 


and a source ~ 

The evaluation of is achieved by the same procedure as for part A, where 

now we use Eq. (|B32|) in (|3^) and (pS]) to solve for 6E[^^ and SE^^^ respectively. 


This computation is somewhat more laborious than that in part A, due to the (3^/3,, 
tensorial dependence and the factor present in this part of the (GM) photon 

propagator. 

To evaluate A, A, and A we need to modify the BBF technique by using ([B5|) 
along with 

(B33) 

to regulate (P32|) . The expressions for the J’s are somewhat more complicated than 
those for 6E^^^ (as expected); but their manipulation and further algebra follow from 
BBF |T^ . The relevant details are in appendix P3.VI1] ; the result for the one potential 
part is 


6E[^^ = 


a 


37rm2 

+ 

(3-p 


7^^ f ^{.p'){ 


17 
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0 m 


(B34) 


+ (- 


m 
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TT 
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y — {3 ■ Vm)ia^^UiQj — mP ■ qia^'^V^Py 




which is good up to order a(Zay, and we have retained only the leading terms as 

/i —0. 

The evaluation of SE^^^ is quite analogous to that for dE^^K The starting point is 
Eq. (P316|) , where (M^) and {M^y are still dehned by (|B2CI|) and ([B28|) respectively. 


We give calculational details in appendix |B.VII|, and quote here only the hnal result: 


SE'P 
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'i'Km? 


7' 




5 ^2 J- A “ M / h 

—vr - h (— I -) ln(- 

12 12 A 2 ^ AW 


(n|VVo|n) 
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(B35) 


We now add 


to (|B35|) to obtain 


5E 


(B) 


a 




r r 11 ^2 47 ,1 ^2 m / 

H “ -717 + (o+“ )ln( 


12 


48 


m 


‘2,Ei^ j 
1 


(n|VVo|n) 


+ ln( 7 —)(n|(M ■ V) Vo|n) + UiUjC"^ + (77 + ^ )C 


'2E, 


(B36) 


y — (3 ■ Vm)ia'^^UiQj + mu ■ qia^^Vfj,(3i, 


,u 


m( 




where we approximated 7 ^ ~ 1 + in order to keep terms only np to order 
As a cross-check on the above resnlt we note that, before expanding 7 ^, the limit 
IdfiPu —*■ Ppiu', yields dE^^'^ —»• —2^'^‘^6E^. This is as expected since according to ( [B32|) , 
G!.? ^ -2«7=G2„. where is the standard (metric) propagator. 

We close this section with a comment on the renormalization procedure. For SE^\ 
the counterterm 6 C was related to mass renormalization. However in this part of the 
calculation we must also account for the renormalization of the THep parameters, 
which show up as functions of: the limiting speed for massive particles, Cg = Tq/Hq, 
and the photon velocity, cl = I/poCq. Charge renormalization is not necessary here 
because the Ward Identity forces a cancellation between the divergences coming from 
the one potential part and many potential part of the self energy. Details of this 
process are shown in appendix [B. VI. 


B.III Vacuum Polarization 

We now need to obtain the vacuum polarization contribution. To the desired approx¬ 
imation, the electrons forming the loop in diagram 1(b) can be considered free. This 
is because Furry’s theorem implies that the next-order correction to this is a diagram 
which contains a loop with 4 vertices, which is expected to be of order a{Zay. In 
that case the result is known to be 

5Ep = j ilj{p')iIlf^'"{q)iG,,^{q)YVf,{q)'ilj{p)d^p'd^p, (B37) 

The evaluation of is identical to the standard (metric) case, since it only involves 
the product of fermion propagators, which are unchanged by the THep action. The 


















differences appear in the renormalization process, where both the charge and the 
Thfe/i parameters must be renormalized, the details of which are shown in appendix 
p.VI| . The result is 




q^qn 


(B38) 


ISvr 

If we substitute Eqs. (p!4D and ( |B38|) in ([B37|) , we obtain after some manipulation 

?2 


6Ep = 


a 


{(n|V'V„|n>( - ^ - |{n|(S' V) V„)|n>} (B39) 


5 


37rm2 L' ' ■ "I ' V 5 
We next proceed to add together the self energy and vacuum polarization contri¬ 
butions to the level shift. 


(B40) 


B.IV The total GM Radiative Correction 

Up to this point we have been able to solve the level shift in terms of 

6En = 5E^^'^ + SEf^ + 5Ep 

where each term has been defined in Eqs. ([B31|) , (|B3(j|) and (|B39|) . 

We note that in 5Es there are terms proportional to 7 , which mix large [ip] and 
small component (y) of tlj. Within the accuracy required we can relate them by 
X = —and so write everything in terms of the large component only. 

Replacing the expression for the external source (^) in ( |B40| ) , we obtain after some 
algebra 

(B41) 


5En = 

“iTim? . 


(1 + e(- + U^))C + + {n\E\n) 


where C and (7*-^ are defined by Eq. 


, and 


E = AT:Za5{x) 


19 , , m 

-h n( 

L30 ^ 


1 58 ^2 /3 2 ^ 2 \, / \ 

2B.'+h-M“45“ + b + 3 “ ' 


Za 


4+^8 


u 

T 


[u ■ n) 


a ■ L 


(B42) 


^^[3{u ■ hf - u'^] 


r-' 
27 T 


14 , m 

-h 21 n(- 

15 ^2E* 


-M - ha -{u xp) — a -{ux h)u ■ p 


We have omitted operators with odd parity (such as m x h ■ a) in ( |B42|) , since their 
expectation values vanish for states of definite parity. 

There is still an implicit dependence on ^ and u in ( [B41| ), which comes from the 
Dirac states (as seen at the end of Sec. III). Note that up to this order all atomic 
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states and energies refered in Eqs. ( |B41|) and (|B25| ) are considered within a non 
relativistic approach. 

In terms of the formal solntion for the Dirac eqnation (p!?!) , we can single ont the 
complete ^ dependence in ( |B41|) , and write 


SE„ = 


a 




{(l + ^d + u^))C°+ + °(n|E|n)°} (B43) 


with 


UiUjE^^ = UiUjC^^ + C' + (°(n|4=o|n)' + h.c.) (B44) 

where C' gronps all the terms in Eq. (|B25|) depending on the pertnrbative states (|n)') 


or energies (E^) as introdnced in Eq. (|^. These pertnrbative states are needed not 
only for the In) state related to the level shift, bnt for all the intermediate states 
introdnced by ([B25|) as well. Eq. (|B43|) is valid np to 0(^)0('U^)0(a(Za)^). 

We can dehne the reference energy E^ as in the metric case by 


En)ln\Er En\ 


-(Za)Dn(^) = E. mmWEr - En) In 


Er—En 


for / = 0 


for / 7^ 0 


(B45) 


where the snbscript 0 has been omitted in the energies and states. This dehnition 
rednces 

-V Q I 0 for / = 0 

^ = \ 4^(Za)Dn(^) for / 7 ^ 0 

which provides an elegant way to write the “Bethe-snm”. The presence of prefered 
frame effects will indnce more “Bethe-snm”-like terms in which, along with the 
contribntion from the pertnrbative states (both ones connted by 6 E^^) will have to 
be evalnated nnmerically for any particnlar state. 

For the Lamb shift states we can use m ■■ 


Ef = 16.640 Ryd Ei'' = 0.9704 Ryd 


2P 


(B47) 


and simplify the last term in Eq. ([B43|) as 






(Za) 


(Za) 


-m 


-m 


19 m rl 58^2 2 ^2 n w 1 

^ “ ^[3 ^ 45“ “^2^3“^ J 


(B48) 
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12 
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107 


- ■= -f cos^ 0 i -h 


eyio 


12 eyio, 


where 6 represents the angle between the atom’s quantization axis and the frame 
velocity u. 
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B.V Many potential part approximations 

In this appendix we justify the following approximations: 



— {^nr^nr^nr) 

(B49) 



(B50) 


~ 0{{Zafa) 

(B51) 


following arguments similar to those presented by BBF . 

We hrst note that, as powers of Za, the orders of magnitude of the different 


terms involved in the expressions in (|B49|) are equivalent to those for the metric case. 
For example, if we look at the source, we see that ~ e(j), where is given by 
Eq. ( 0 ) and 0 is the ordinary Coulomb potential, and so the relative order between 
the nonmetric and metric case is the same. Furthermore, as discussed at the end 
of Sec. Ill, the states |n) and |n)° also have the same order of magnitude, as do 
the quantities En and E^. Discrepancies that could be expected from the photon 
propagator, particularly from the part proportional to (3^(3'' (in contrast to the 77^1/ 
dependence for the standard case), are not important as long as the transversality 
condition is satished for the M operators, since this condition relates the differing 
components with the appropriate orders of magnitude. Finally, unlike the photon 
propagator, the bound propagator retains the same form as in the standard case, with 
differences arising only from the expression for the external source. As a consequence 
its further simplification is analogous to the metric (BBF) case. 

Let us look at the many potential part. From (^) we get 


{MKlM) = 

X K^{En - ko;p'- s'- k,p +s-k) 
Ml{p + s- k,p)iJn{'P}G^'"{k) 


(B52) 


X 


for the generic structure of the terms on the left hand sides of (|B49|) - (|B51|) , where 
the constant factors and integrations over pi and Sj have been omitted. The nonrel- 
ativistic and relativistic regions are dehned according to \k\ ~ {ZaYm « m and 
\k\ > m,respectively. In considering the relevant orders of magnitude in each of the 
expressions (|B49|) - ([B51|) that follow from ([B52|) , we note that, to lowest order in Za, 
the relevant contribution from comes when fco ~ 1^1, and that we can employ 
the nonrelativistic expressions for the 'ipn, making use of the approximations given by 


(H) 
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Turning now to the relation ([B49|) , we can prove it by showing that the contribu¬ 
tion of relativistic states for is of a higher order of magnitude than for . We 
can see from ( pi2|) and (|B13|) that differs from by a factor (leaving aside the 
temporal component) {ff which in the relativistic region [ko ~ m) is of order 

Za. Therefore the contribution of in that domain will be of at least one order 
higher than that of . Since the latter is already of the desired order (assuming 
the validity of ( |B50| ) ) we can neglect the contribution of the relativistic states for 
and consider it, along with the bound propagator, in its nonrelativistic limit. 

To prove the relation (P50|) we evaluate the error due to the neglect of the electro¬ 
magnetic potential in the intermediate states. We imagine that one extra potential 
{y) acts between and . This introduces an extra factor of order 




drr 


-r) 


— K — m 


(TV 


2kom 


^ ~ {Za) 


(B53) 


which is negligible within the accuracy required. We have then shown that, in the 
evaluation of the intermediate states may indeed be regarded as free. 

The relation (|B51|) follows from arguments similar to those used to justify 


Since in the relativistic region is one order higher than , the cross term in 

that region will also be one order higher than and so is negligible. On the 

other hand in the nonrelativistic region will be dominant (note the factor /cq in 
its denominator) over . That is 




f? — p 


Za 


(B54) 


and so the product of these terms will be negligible in comparison with (M^). Hence 
the cross terms yield results that are at least one order higher than the desired order, 
and so they do not need to be included. 


B.VI Renormalization 

Just as in the standard (metric) case, we need to renormalize the various parameters 
of the theory in order to get rid of the divergences. In the standard case, those 
parameters are the mass and charge, although the latter only needs to be renormalized 
for the vacuum polarization contribution. The self energy part has no need for such 
a renormalization, since the divergences coming from the one potential and many 
potential parts cancel each other. In the nonmetric case, we have also to include the 
renormalization of the THefi parameters, which show up as functions of Cq = Tq/Hq 
and c* = l//ioeo- 
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In part A of the calculation, renormalization is identical to the standard case. 
The counterterm 6 C is just related to mass renormalization. In part B, we need to 
consider additional counterms, since 5C should also account for the renormalization 
of the T Hefi parameters. 

In units where cq = 1 (c* = c), EEP-violating corrections only appear in the 
electromagnetic sector of the action (as terms proportional to ^). However we could 
choose more generally cq 7^ 1, for which the particle sector of the Lagrangian density 
is of the form 

Cd = Y - m)Y + ioYiVo - Ao)7V (B55) 

with ,^0 = 1“ Cq or in the moving frame (after using ([T0|) ) is 

C'^ = (B56) 

+ -p-P-V) 


up to a constant. 

From (|B56|) we see that quantum corrections of the form 




( 2 ), 


(B57) 


can still be expected. Note that gauge invariance will guarantee = S^q. 

Hence, in order to renormalize the mass and the THe^ parameters, we have to include 
counterterms of the form 


5C = 6 m + 6^0 0(/3 ■ p — /3 ■ V) (B58) 

where 6 m and are chosen such that 6 Bs gives zero contribution as the source is 
turned off. This condition forces /a = 0 when acting on free spinors. 

Finally, for the vacuum polarization contribution the charge has to be renormalized 
along with the THep parameters. Charge renormalization is identical to the standard 
case. For the THefi parameters the procedure is equivalent to the self energy part, 
where now, given the form of the electromagnetic action (see Eq. (0)), we expect 
quantum fluctuations of the form 

6Cem = 6^A>^{{k^ - iP ■ kY)p^, - P^P.k‘^}A- (B59) 

to occur. Hence a counter term of that form it is needed to renormalize the THep parameters, 
or equivalently .^ = 1 — Ho/Topoeo. 
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B.VII Calculational Details of Type B Contributions 


We present here further details underlying the computation leading to Eqs. 
and (|B35|), which are referred as the type-B contributions to the self energy. In this 


part the photon propagator to be considered is given by ( |B32| ), where the hrst and 
second terms have respectively a tensor dependence like and and need to be 
regularized according to ([B5|) and ([B33|). We show the relevant details involving the 


hrst term of the propagator only, since the remainder can be computed in a similar 
way. 

We begin then with the one potential part by simplifying Ji. After replacing 
in (0. we get 


h = 




2p' ■ (3— ^ ^ 2p ■ f3— ^ 0 d'^k 


f 


-dL + 


(B60) 


J k? — 2p' ■ k k3 — 2p ■ k {k^ — Ly 

where from now on ■ ■ • stands for the contributions coming from the second term of 


If we use 


1 ^ef'dx 


^(1 — z)dz 


abc^ Jo Jo [{ax + b{l — x)){l — z) + cz] 
we can rewrite Eq. (P360|) as 


(B61) 


h= - Ap-yp'-yfJo + 2 p-( 3 y^^fjy 


where 

with 


+ 2p’-i3fY fi 

dLd'^k 


+ ■■■ 




[{k-pyi - z)f - Al] 


-{l]k^]k^ky] (B62) 


h = 


Px = xp' + (1 - x)p Al = pI{ 1 - zY + Lz 
After evaluating (P62|) , we can express 

[A ■ pO ■ p'(ln 4 - 2) + f p^(| - In ^)] 

J Pi^ ^2 /i^ 8 2 Pi i 

x^{fj-pr Y+fj-p' yr)-{fj-pr Y+fj-p' yr){i-x)^ 

4 


(B63) 


+ yx[fJ-V{(J-p +(J-p')-^Px-V) 

+ ■ YPx-V -^V ■ I3pl{l + ^ - In ^ 

— (1 — x)P ■ pV ■ p — xp ■ (JV ■ I + ■ ■ ■ 


(B64) 
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The evaluation of the remaining J’s is analogous, and so 
^2 ^2 1 „2 

h = -T/3'P^7"5(ln^ + |)-ij<;3V5(ln^-h + "- + '5C' (B66) 




(B65) 






From appendix |B.VI| , we know 6 C = 6 m + 6 ^o 13{P ■ p — P -V), where in this case 


6 m = — 7 ^^(ln ^ “ “) + 
8 m^ 


1 2^/, A2 5n 

= p + 5 ) + 


(B67) 


Since here = rj'^^Vo, we can rewrite after some manipulation 

Ii + l 2 + h = + K, + K,) + --- (B68) 

where 

^ , , \ (iP-P? 


Ki 

K 2 

Ks 


dx 


pI 
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(my + P-p 


(..y) + i 




P ■ pP ■ p^ + ]^P ■ pP ■ 


1 Tfl ^ 1 Tt^ 

-j \n(—)S(p -p')-^--l3-V 0H^) 

4 p^ 2m 4 m^ 


1 

2 


dx 


pI 


-^xp ■ p' p' p y + {1 — x)p ■ p y p p 
+ r ■ p 1Y + ,3 • Pi h - 2/3 ■ 17(;3 ■ p + ,3 ■ p') + 1,3 ■ p / >7 ^ 

} 


{i-x)p-p' isr'^- pI pf-2p-vp+rp— 

1 m 


We want a result good to a^Zap, and so we can simplify the above expressions 
by using the assigned order given by (p3|) , from which we can relate 


q = p' — pr^ Zam 

p'2 _ p2 ^ . p^2 _ p2 ^ p2_^2 ^ [Zapm^ 


and then reduce to 


K. ~ 


JL 

2 m^ 




P \ ^ J2t^ I o2\ 


H^)-lq\^ + PP + {P-q? 

mb! 


(B69) 


(B70) 


where antisymmetric terms under p' ^ p vanish. 

To simplify K 2 we follow BBF and use 

{p-mf p = y{p-m) p = 2p-py-2my p- y p y 
~ 2{y-p-my) p-4p y-v^ 


(B71) 
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where we have assumed the operator is acting on Dirac spinors of momentum p and 
omitted the integration coming from 

i){p){]^-m) =6{qo) J 'il){p') f{q)d^p' (B72) 


Note that y ^ ~ since the square of the potential (after factoring out 

the spinors and integration variables) is already of the desired order (Za^ (see (p^) ) 
and so ~ 7 o — 1. 

The hnal result is 


Following a similar approach we reduce 

Ks ~ -^(V-p-mr) + — + /3-V^-^ ^V-p (B74) 
2 m m m 


^■<1 .y_P-y 


2 m 2m 

We can make further simplihcations by using 


3^ 3"^ 


'ip{p')B{p', p)'ip{p)d^p' d^p = 0, 


(B75) 


provided A{p', p)A = —B{p,p'), where B represents any operator as a function of 
p' and p, as for example, 3' Q Y■ Note that we are interested only in the real part of 
the level shift. 

Putting everything together, we obtain after some manipulation 


= 


a 


Tim 


Ai! Y{p'){f{3■ A 1 + Ih^ 


m 


fq^ ^ + - + -ln(^) 
^^ 16 8 6 


f + 3- V^)ia^^Uiqj - ^3 ■ qia^^''Vf,3u 


(B76) 


+ - ]^)i(^^''V^qy]3{p)d^p'd^p - + ■ ■ ■ 


Note again that this represents the calculation involving only the hrst term of Eq. 


Now to evaluate the many potential part contribution we need to solve Eq. ([B16|) , 
with {m3 and {M^3 given by Eqs. (|B2CI|) and ([B28|) respectively. 

So, after substituting ( [B32| ) in ( |B20|) 


{m3 


a 


Ati^i 




\3-M3 d^k 

-fjfi k^-En- Er ^ ' 


(B77) 
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with 


= (n|i?^|r), 

Using the transversality condition, we relate 

Kq Kq 

which reduces the integral on the angles of k to 

P- 

dVt\(3 ■ + \u ■ M 

OfCn 


(B78) 


We evaluate the remaining and \k\ integrations in ([B77|) , by using ([B21|) , ([B23|) 
along with the analogous relations 


u- M. = 


-1 

m/co 
1 


{Er - En){n\u-^r) 


^\{n\u-^r)\^{Er-En) = -{n\u-Wo\n) 


to hnally obtain 


a 


(M^) = —I (7 + 

7rm n 


1 A 1 

_ , UiU. 

6 2 ^ 


UiUjC^^ + 


2 1, / N 

- + -ln(-^) 
9 6 ^2EJ, 


(n|VVo|n) (B79) 


fi 


7 + 7in(^) 

2 2 ^2EJ, 


(n|(M-V)Vo|n)} + --- 


where we have kept only the leading terms as /i —> 0 and neglected the imaginary 
part. 

The computation of (M^^) is straightforward. Here we need to replace (|B32|) in 
(|B28|) , and use Vk = VaoV^- Further simplihcations follow from BBF and the assigned 
order of magnitude given before. The hnal result is 


(M^^) = + ■ ■ ■ 


a 




*^0 


TT 


4m 


(B80) 


Adding together 


, and ([B80|) will give us then the final expression 


for the self energy contribution for this part of the calculation. Note that the above 
results can be verified by taking the limit > 7 ^,^, which reduces 

and therefore the former expressions should reduce up to a constant, to the metric 
case. 
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C Virtual non metric anomaly 


In the THefi formalism, gravity interacts with matter through the T and H functions, 
which are assumed locally constant within atomic scales. A priori they do not need to 
be the same for differents type of matter (like baryons and leptons), or furthermore for 
matter and antimatter. In this context for example, a non metric anomaly related to 
electron/ positron difference will modihed the Lagrangian density related to fermions 
by 

Cd = Y - v[i)i} + - Ao) 7V^ (Cl) 


where = 1 — c_/c+ and Cip = with — and + labeling electrons and 

positrons respectively. After using ([T0|), we can refer (|CI|) to the moving frame as 


C^D = Y{i>- Y - m)ijj + ^+7V"^(/5 -p- f3 -V) (C2) 


The imposed broken symmetry between particle and antiparticle changes the 
fermion propagator (in the positron case) to (up to 0(^+)): 


Sp = {'^ — m) ^ + ^7^ ■ p(^ — m) ^ (C3) 


where the hrst term represents the unchanged electron propagator S^- 
The positron-electron pairs produced in the electric held of the atomic nucleus, 
are seen in the Lamb shift transition via the vacuum polarization contribution given 
by ( |B37| ), where in this case: 

av“(q) = ^(-1)^’’- / dWiS^ip + qh'-^s+ip) (C4) 


After using eq. ( PS]) along with standar technics p7 
metric part of ( |C4| ) is up to O(g^) 


we obtain that the non 


a 


27r 




{4,)=} + ... (C5) 


where ■ ■ ■ accounts for the gauge dependent terms which give no contribution to ([B37|) . 
Eq ( |C5| ) also comes after proper regularization and renormalization processes, which 
follow from previous sections. 

In this EEP violating context, the radiative corrections related to atomic energy 
levels are modihed by (up to 0 {a{ZaY0{uY) 


SE+ = SE+ = + {n|(S. V)=V'„)|n)} (C6) 


43 










where we have replaced in ([B37|) and simplihed afterwards. By taking the Lamb 
atomic states, we hnally obtain 


777 

^Et = (C7) 
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